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Abstract. In cool stars that oscillate like the Sun, non-radial modes become mixed 
QQ ' as the stars evolve. The mixing is caused by the coupling between g-modes in the stellar 

core and p-modes in the envelope, which results in distinctly different and more com- 
plex frequency spectra for subgiants and red giants than seen in main sequence stars. 
, Using a new version of the 'scaled' echelle diagram, I illustrate how the frequencies of 

Pm ! non-radial modes evolve during the evolution from the main sequence to the red giant 

jy^ • branch, and I show how they depend on stellar mass and metallicity. Then, with focus 

(-H ' on the dipole (1=1) modes, which show the strongest effects from mixing, I present a 

O I toy model to fit, and hence identify, those modes in a large series red giant models. 

o 

C^ . 1. Introduction 

^ . Stellar models predict that non-radial oscillation modes in stars lik e the Sun w ill be- 

K»- ' come of a mixed character as the stars evolve off the main sequence dOsakill 19751) . The 



mixing, which is caused by the coupling of two modes (a g-mode in the stellar core 
and a p-mode in the envelope) will result in shifts of the mode frequencies - known as 
mode bumping - as the two modes undergo so-called avoided crossings. As a result, 
mixed modes do not follow the usual asymptotic rel ation of regular spaced high radial 
order p-mode frequencies, v , as expressed by (e.g. IVandakurov.1968. : Tassoul. 19801 : 



Christensen-Dalsgaardll201 ll) 



Vni = Avin + 1/2 + e)- Svqi (1) 

K^ ^ where the small frequency separation 

H ! 9 Ay 

. .C^; Svoi^il+l/2f—. (2) 

Vnl 

Here n is the radial order, / is the degree, e is a dimensionless offset, and Av is the sepa- 
ration between over tone modes of the same degree. When visualized in the echelle dia- 
gram, where we plot mode frequency versus the frequency modulo the large frequency 
separation. Ay, mixed modes can show strong departures from the almost vertical ridges 
of the non-mixed modes, which I will discuss next. 

2. The scaled echelle diagram 

A useful diagram to compare r node frequencies of many st ars (or models) is the scaled 



echelle diagram introduced by lBedding & KjeldsenI (l201Clh . Their Method 2 sought to 
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Figure 1. Asteroseismic Hertzsprung-Russell diagram with v^ax as luminosity 
proxy. Tracks are for Z - 0.017 and their masses in solar units are shown at the 
ZAMS. 



make the ridges of the radial modes coincide for different, though still similar, stars 
by fine tuning Av; in essence removing any differences in e by allowing the ridges to 
depart slightly from vertical. This provides a clear picture of how the modes of different 
degree are positioned relative to one another in a comparative analysis. Here I will 
use a slightly modified version of this method to show the evolution of a large series 
of models ranging from the main sequence to the red giant branch (RGB). For each 
model, Av is derived from the four radial modes clo sest to the frequency of max imum 
power, Vmax (scaled from the solar value according to lKjeldsen & BeddingI (Il995b ). and 
the variation in e from one model to the next is removed by shifting the frequencies 
to a common value of e = 1.2; hence maintaining strictly vertical ridges of the radial 
modes. This could therefore be appropriately called a shifted echelle diagram. 

The stellar models were derived with ASTEC dChristensen-Dalsgaardl 2008al) and 



adiabatic frequencies were calculated using the ADIPLS code (IChristensen-Dalsgaard 



2008bl) . The models range in mass from O.SMq, to 1.6Mq in steps of 0.2Mq, and in Z 
from 0.01 1 to 0.028 in steps of 0. 1 dex in log Z/X (with X = 0.7 fixed). The tracks were 
evolved from the zero-age-main-sequence (ZAMS) towards, bu t not quite reaching , the 
tip of the RGB. Details on the model physics, can be found in IStello et al.l (12009 ') and 
references therein. The tracks for Z = 0.017 (approximately solar metallicity) are 
shown in Figure [TJ where I have used Vmax as the proxy for luminosity to make it more 
useful in comparison with the following figures. 

In Figure [2l I show an echelle diagram for a series of models along the one solar 
mass track, where four orders are plotted for each model. The symbol sizes are scaled 
according to the reciprocal of the square root of the mode inertia, normalized to the 
radial modes, but set to zero below a certain threshold for clarity. Several features are 
apparent in this diagram. 

1. We see an increasing tilt of the non-radial ridges with increasing /. To illustrate 
this, I have connected modes of the same degree by lines for the ZAMS model. 
The tilt can be explained by the last term in the asymptotic relation, dvoi (Eq. |2ll. 
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Figure 2. Echelle diagram for series of models along a I.OMq track (Z = 0.017, 
[Fe/H] ^ 0.0; see Figure [TJ- For each model the four modes closest to v^ax are 
plotted after they were shifted to a common value of e = 1.2. The dotted line shows 
the halfway point between adjacent radial modes. 
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where a higher I implies a higher frequency dependence, and hence a stronger 
tilt. 

The position of the non-radial ridges, and hence the small separations, changes 
quite significantly with evolution. After moderate changes during the main se- 
quence, there seems to be a slight pause during the subgiant phase where the 
stellar density is not changing by much, before more dramatic changes occur 
during the ascent of the RGB. That the change is larger for increasing mode de- 
gree is explained by the same asymptotic term as for the above tilt. Interestingly, 
this evolution of the ridge positions shows s ome similarities to the variation seen 
in e as a function of Vmax (IWhite et al.ll201ll) . 



At a certain point during the evolution, corresponding roughly to the transition 
between the main sequence and the subgiant phase, mixed modes start to appear, 
most clearly for Z = 1 and 1 = 2, leading to modes spreading across the entire 
echelle diagram (Z = 1 in particular). 

After the subgiant phase the mixed modes become more confined as the models 
evolve up the RGB. Ultimately, neai" the tip of the RGB, the coupling between p- 
and g-modes become so weak (efficient mode trapping) that we expe ct to see only 
one m ode per degree per radial order, as in the main sequence phase (iDupret et al.l 
2009h . 



5. As first observed by lBedding et al.l feOlO) and shown for stellar models by Mont alban et al 



(i2010,) . the small separation, 6v\o, between the the Z = 1 modes and the point 
halfway between consecutive radial modes (dotted line) becomes negative for 
the red giant models. 

Note that the less dense 1 = 5 ridge arise simply because these modes were not 
derived for all models along the evolution track at the subgiant and red giant phases. 

2.1. Changing mass and metallicity 

To illustrate how the different features seen in the echelle diagram depend on stellar 
mass and metallicity, I show examples for Z ^ 3 of models with a range of masses and 
metallicities (Figures |3l IH). Due to the high computational requirements for deriving 
mode frequencies of evolved stellar models there is a lack of frequencies in the model 
grid, which appears near the base of the RGB (300//Hz < v^ax ^ SOOjuHz). Fortunately, 
the highly evolved red giant models at even lower y^ax enables 'interpolation by eye' 
to bridge the gap. 

It is evident from Figure |3] that the ridge positions change significantly as a func- 
tion of mass, with lower masses showing smaller separations between the odd ridges 
iSvo2) and even ridges (Svu), respectively. The higher mass models, which exhibit the 
smallest density changes during the subgiant phase, show more static ridge positions 
during this phase in agreement with Eq.[2]and the observation made in Figure [2] (point 
2). We further see a hint that the Z = 2 modes spread out more for the least massive 
models, indicative of stronger mode coupling. 

If we turn our attention to the effect from metallicity (Figure ^, we see that an 
increase in Z results in a decrease in the separations between the odd ridges and even 
ridges, respectively. This is in qualitative agreement with the effect we saw on mass, 
namely that the cooler the model the smaller the separations 5vq2 and <5vi3 become. 
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Figure 3. Echelle diagram for series of models along the A: O.8M0, B; 1 .2M0, C: 
I.4M0, and D: I.6M0 tracks (Z = 0.017, see Figure[T]). The y-axes are rescaled for 
the lower panels for improved visibility. 
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Figure 4. Echelle diagram for series of models along 1 .OMq tracks of A: Z = 
0.011 ([Fe/H] ^ -0.19), B: Z = 0.014 ([Fe/H] ^ -0.08), C: Z = 0.022 ([Fe/H] ^ 
0.11), and D: Z = 0.028 ([Fe/H] ^ 0.22). 
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Figure 5. Freque ncy spectra of an RGB star and a red clump star (taken from 
iBedding et al.ll201 U) . Mode degree is indicated. 



2.2. Mixed dipole modes in red giants 

In the following I will focus on the lower part of Figure |2] where we see increased 
complexity in the 1=1 ridge due to mode mixing. In particular I will address how to 
identify the 1=1 modes in frequency spectra of red giants. 

The complicated frequency spectra of Z = 1 modes in subgiants can i n many cases 

be su ccessfully fitted, resulting in correct identification of the mode degree (IDeheuvels & Michel 
20101 . Benomar et al. in prep.). However, can something similar be done for red giant 
stars? To motivate this question, I show the frequency spectra of two red giants in Fig- 
ure \5\ The weaker mode coupling in red giants imply we do not see the 1=1 modes 
appear everywhere, as in subgiants (Figure O, but rather within well confined clusters, 
one for each radial order (centered at the frequency of the 'unperturbed' p-mode in the 
hypothetical case where there would be no coupling with the g-modes in the core). On 
the one hand, this makes the spectra more simple as we avoid a very dense spectrum of 
1=1 modes that overlaps with all the other mod es. On the other hand, the slightly com - 
plicated frequency separations of these modes (iBeck et alJl201ll : lBedding et alj|201ll) 
makes it difficult to estimate how many modes are in fact 'missing in the gaps' between 
each cluster of modes, and hence hinders a straightforward determination of the fre- 
quencies of the undetected modes. In practise this means that it can become difficult to 
make judgment on whether a peak in the observed frequency spectrum belongs to the 
sequences of / = 1 modes or whether it is due to a mode of a different degree, or simply 
noise. 

Fortunately, for stellar models we have frequencies for all modes available, even 
those in the gaps that are not detectable in reality. This can help us develop ways to 
analyse the complex frequency spectra observed in red giants. As a first step towards 
fitting observed 1=1 modes I will show some preliminary results for what one ob- 
tains when fitting a simple toy model to the stellar model frequencies. The aim of this 
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Figure 6. Period spacing versus frequency for dipole modes (diamonds) of Model 
A, which has T^ff = 4745 K and L/Lq = 56 (marked in Figure |7j5). Toy model fit 
is shown by thick solid blue line. Resulting fitting parameters, APg, FWHM, and 
Depth are indicated (see text for details). 



exercises is to identify relations between on the one hand the parameters that describe 
the toy model and one the other hand the quantities that we are able to detect in the 
observed frequency spectra - ultimately enabling us to make informed decisions and 
on how to interpret the observations. 

Figure [6] shows the period spacing of the 1=1 modes from a model located on 
the RGB. The period spacin g is approximately constant as expected from asymptotic 
theory of pure g-modes (e.g. 'Christensen -Dalsgaardll201 iL and references therein). The 
departure from the constant spacing, seen as regular dips, are caused by the coupling 
between the g- and p-modes, which shifts the frequencies. The dips are separated by 
Av- the separation between the unperturbed p-modes. The most weakly coupled modes 
are almost pure g-modes, hence showing near equal period spacings. However, being 
pure g-mode means they are basically contained in the core, and hence not observable. 
The more strongly coupled modes, which show the largest departure from equal period 
spacing, are more p-mode like, and have therefore significant amplitudes at the surface 
and are the ones we detect as small clusters. 

I have simultaneously fitted a series of Lorentzian profiles, one for each dip to a se- 
ries of stellar model frequencies of dipole modes. The free parameters in the fit are AP^ 
(the true g-mode spacing unperturbed by mixing with p-modes), FWHM (the width of 
the Lorentzian), 6{FWllM)/6f (the change in the width as a function of frequency). 
Depth (the depth of the dip), 6(Depth)/6f (the change in the depth as a function of 
frequency), and the frequency of the centre of each Lorentzian profile, assumed to be 
strictly equally spaced. Figure |7^ shows a 2.4M0 track evolving from the base of the 
RGB to the base of the asymptotic giant branch (like a transposed Hertzsprung Rus- 
sell diagram), and Figures |7j5-d show the result of fitting the toy model to the stellar 
models along that track. An example of a fit is shown by the blue curve in Figure |6] 
Interestingly, there seems to be a strong relation between APg and Depth, but Depth is 
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Figure 7. (a) Evolution track of 2.4M0 track, (b)-(d) Results of fitting the toy 
model (Figure |6] thick blue curve) to each model along the track. Plotted values of 
FWHM and Depth are evaluated at Vmax- 
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not simply jAPg as suggested by iBedding et all (|201 ih . Although FWHM and Depth 



are both related to the coupling strength between the g- and p-modes, only Depth sep- 
arates the RGB from the red clump, while FWHM is largely a function of Teff. When 
repeated for stellar models of different masses, the relations between the toy model 
parameters and the seismic observables change. Quantifying these relations will help 
create a formalism to fit mixed mode frequencies from observations. 



3. Summary and future work 

I have shown that a new version of the 'scaled' echelle diagram can be useful to com- 
pare an ensemble of stellar models, and in particular to show the evolution of pulsation 
frequencies from the main sequence to the RGB as well as how the frequencies depend 
on mass and metallicity. For red giants we see that the 1=1 modes become more 
confined as small clusters of modes. While this in some sense makes the frequency 
spectra simpler, it also brings new challenges by making it harder to determine the fre- 
quencies of all the modes between each cluster. I briefly demonstrated that fitting a 
simple toy model to stellar model frequencies can aid in establishing relations between 
seismic observables and the toy model parameters. When fully developed this could 
lead to a robust way to fit the observed 1=1 frequencies and hence to determine the 
characteristics of their complex frequency separations. 
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